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1  Abstract 

2  The  factors  contributing  to  the  timing  and  intensity  of  phytoplankton  blooms  ("red  tides") 

3  in  very  shallow  water  are  poorly  understood.  Here,  the  sudden  and  brief  (3  h) 

4  appearance  at  the  surface  of  an  alongshore-parallel  band  of  red  tide  is  described.  The 

5  bloom,  near  Huntington  Beach  CA,  was  dominated  by  the  locally  common  red-tide 

6  forming  dinoflagellate  Lingulodinium  polyedrum  (F.  Stein).  Surface  chlorophyll-a  (Chi  a) 
1  and  temperature  (7)  were  mapped  with  a  novel  GPS-tracked  jetski.  The  surface  red-tide 

8  band,  with  maximum  Chi  a  ~7  pgL'\  was  visible  for  more  than  3  h,  was  located  500  m 

9  from  shore  (~8  m  water  depth),  was  -200  m  wide  and  extended  -1  km  alongshore.  The 

10  subsurface  bloom,  with  maximum  Chi  a  >30  pgL"',  developed  farther  offshore,  near  the 

11  17  °C  isotherm  over  4-7  days  prior  to  the  surface  appearance.  A  few  hours  before  the 

12  surface  band  appeared,  an  intense  Chi  a  patch  was  observed  in  13  m  total  depth  in  the 

13  mid-water  column  trough  of  a  shoreward-propagating  supertidal  internal  wave.  The  Chi  a 

14  patch  intensification  was  consistent  with  dinoflagellate  vertical  swimming  interacting 

15  with  the  wave-driven  circulation.  In  shallow  water  (<12  m  depth),  vertical  mixing, 

16  attributed  to  internal  wave  breaking,  brought  the  Chi  a  patch  to  the  surface.  Warm  water 

17  in  very  shallow  water  (depth  <  5  m)  apparently  blocked  the  Chi  a  patch  from  entering  the 

18  surfzone.  These  observations  highlight  the  role  of  small-scale  physical  and  biological 

19  processes  in  determining  the  spatial  and  temporal  characteristics  of  nearshore  red  tide 

20  blooms. 
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1  Introduction 

2  Intense  phytoplankton  blooms  (“red  tides”)  and  harmful  algal  blooms  (HABs)  in 

3  the  nearshore  (<15  m  depth),  incur  losses  of  millions  of  dollars  annually  to  aquaculture 

4  (Alonso-Rodriguez  and  Paez-Osuna,  2003),  threaten  marine  mammal  and  human  health 

5  (Anderson,  1997),  and  inhibit  beach  recreation  (Backer  et  al.  2003).  Nearshore  red  tides 

6  are  difficult  to  predict  and  monitor  (Anderson,  1997),  in  part,  because  phytoplankton  are 

7  transported,  mixed  and  concentrated  by  a  range  of  physical  processes,  including 

8  upwelling  (e.g..  Smith  et  al.  1983),  winds  (e.g.,  Keafer  et  al.  2005),  alongshore  plumes 

9  (e.g.,  Franks  and  Anderson,  1992),  internal  tidal  bores  (e.g.,  Leichter  et  al.  1998),  internal 

10  waves  (Lenert-Cody  and  Franks,  2002)  and  tidal  currents  and  stirring  (e.g.,  Cloern  and 

1 1  Dufford,  2005).  In  the  surfzone  (within  the  shallow  region  of  surface  wave-breaking), 

12  currents  (Campbell  and  Bate,  1988)  and  species-specific  behavior  (Talbot  and  Bate, 

13  1986)  appear  to  control  the  distributions  of  several  sand-dwelling  diatom  species.  This 

14  broad  range  of  mechanisms  can  create  surface  patches  of  red  tide  organisms  that  may  be 

15  much  smaller  (meters)  and  shorter  in  duration  (hours,  e.g.,  Franks,  1997)  than  the  scales 

16  of  the  growth  and  decline  of  the  bloom  itself  (kilometers  and  weeks,  e.g..  Welters  et  al. 

17  2003). 

18  The  impact  of  a  red  tide  on  the  nearshore  depends  on  the  local  bloom  intensity 

19  and  duration.  Mussels  briefly  exposed  to  high  doses  of  HAB  toxins  may  surpass 

20  quarantine  limits  within  one  hour  (Bricelj  et  al.  1990),  indicating  that  sudden  changes  in 

21  red  tide  concentration  and  distribution  may  carry  important  consequences  for  mussel 

22  health.  “Cryptic"  blooms  may  develop  as  thin  subsurface  layers  (McManus  et  al.  2008), 

23  with  no  (or  limited)  surface  expression  and  therefore  invisible  to  satellites  and  visual 

24  inspection.  In  situ  sampling,  including  subsurface  waters,  at  scales  fine  enough  to  resolve 
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transient  and  subsurface  Chlorophyll  a  (Chi  a)  or  HAB  layers  is  critical  for  effective 
monitoring. 

Many  marine  phytoplankton  are  motile,  and  capable  of  maintaining  position  in  the 
water  column  position  despite  negatively  buoyant  cell  bodies  (Cox,  1980).  A 
combination  of  motility  and  convergent  advection  may  lead  to  accumulation  of 
phytoplankton  at  fronts  (Franks,  1992)  or  in  internal  wave  troughs  (Lennert-Cody  and 
Franks,  1999).  Swimming  has  been  implicated  in  the  accumulation  of  larvae  at  internal 
tidal  fronts  (e.g.,  Clancy  and  Epifanio,  1989)  and  shoreward-propagating  bores  (e.g., 
Pineda,  1999).  However,  the  relative  importance  of  swimming  and  physical  mechanisms 
in  concentrating  and  transporting  phytoplankton  to  the  nearshore  is  poorly  understood. 

In  this  paper,  biological  and  physical  observations  of  a  nearshore  red  tide  are 
presented  from  Fall  2006  as  part  of  the  Huntington  Beach  experiment  (HB06).  This  site 
was  selected  in  part  because  of  chronic  beach  water  quality  problems  (Boehm  et  al. 
2002a).  Previous  studies  include  investigations  of  the  internal-tide  driven  cross-shore 
transport  of  sewage  effluent  (Boehm  et  al.  2002b),  subthermocline  water,  and  suspended 
or  dissolved  material  (Noble  et  al.  2009).  Subtidal  internal  bores  (with  durations  of  2-5 
hrs)  may  also  be  responsible  for  cross-shore  sediment  transport  (Noble  and  Xu,  2003). 
This  is  the  first  study  to  explore  the  high  spatial  and  temporally  resolved  patterns  and 
transport  of  a  red  tide  at  this  site.  Santoro  et  al.  (2010)  use  Sea-viewing  Wide  Field-of- 
view  Sensor  (SeaWiFS)  to  contrast  red  tide  blooms  in  2005  (an  anomalous  year  with  an 
intense  6-month  bloom)  and  2006  (with  comparatively  low  SeaWIFS  Chi  a)  with  a  focus 
on  Huntington  Beach.  However,  satellite  measurements  cannot  detect  subsurface,  short- 
time,  or  small-scale  blooms. 
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1  Here,  the  sudden  surface  appearance  of  a  subsurface  “cryptic  bloom”  of  the 

2  dinoflagellate  Lingulodinium  polyedrum  is  described.  The  surface  expression  of  the 

3  bloom,  visible  from  the  beach  as  a  brown  streak  ~1  km  from  shore  and  extending  >1  km 

4  alongshore  (Fig.  1)  developed  over  the  course  of  a  few  hours  on  12  Oct  2006.  Intensive 

5  and  unique  biophysical  observations  with  high  temporal  and  spatial  resolution,  spanning 

6  4  km  cross-shore,  from  the  surfzone  to  25  m  total  depth  (Fig.  2),  indicate  that  the  bloom 

7  developed  in  the  nearshore,  within  a  subsurface  layer  over  about  4  to  7  days.  An  intense 

8  subsurface  Chi  a  patch  formed  within  a  shoreward-propagating  internal  wave  trough, 

9  consistent  with  models  of  depth-keeping  swimming  behavior  in  the  internal  wave- 

10  induced  circulation.  Internal  wave  breaking  in  shallow  water  (depth  <  12  m)  is  implicated 

11  in  the  vertical  spreading  of  the  intense  patch,  forming  the  visible  band  of  surface  Chi  a 

12  just  seaward  of  the  surfzone,  without  entering  the  surfzone. 

13 

14  Methods 

15  The  nearshore  HB06  field  experiment  was  conducted  at  Huntington  Beach,  CA 

16  from  15  Sept,  to  17  Oct  2006,  within  the  framework  of  the  larger  scale  Southern 

17  California  Coastal  Ocean  Observing  System  (SCCOOS)  and  United  States  Geological 

18  Survey  (USGS)  programs.  Here  the  focus  is  on  a  few  days  of  HB06  observations, 

19  spanning  0.5  km  alongshore  (y)  and  4  km  offshore  (x)  to  25  m  depth  (Fig.  2).  The  mean 

20  (tidally-averaged)  water  depth  is  denoted  as  H,  and  the  vertical  coordinate  z  is  positive 

21  upward,  with  z  =  0  m  at  the  tidally  averaged  surface. 

22 

23 

24 
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Jetski  surface  maps  and  CTD+Chl  a  Casts 


2  Quasi- synoptic  maps  of  near- surface  Chi  a,  turbidity  and  temperature  (7)  were 

3  constructed  at  45  min  intervals  between  1200  and  1500  h  from  data  acquired  by  a  flow- 

4  through  sampling  system  mounted  on  a  GPS-tracked  jetski  (Clark  et  al.  2009).  The  jetski 

5  was  driven  on  1  km  long  cross-shore  transects  between  the  shoreline  and  ~12  m  depth 

6  with  8  transects  distributed  over  600  m  alongshore  (dashed  box  in  Fig.  2a).  CTD  (Seabird 

7  19)  and  flow-through  Chi  a  fluorometer  (WET  Labs  WETStar)  casts  were  made  at 

8  stations  on  a  single  cross-shore  transect.  Cross-shore  transects  were  made  between  H  = 

9  40  m  to  //  =  5  m  (spanning  ~6  km)  from  0900  to  1 130  h  once  per  day  on  various  days 

10  throughout  HB06.  In  addition,  on  12  Oct  at  45  min  intervals  between  //=13mto//=4 

11  m  from  1200  to  1500  h.  Vertical  profile  data  were  averaged  into  1  m  vertical  bins.  Jetski- 

12  derived  surface  transects  {x,y)  and  CTD-i-  Chi  a  -derived  cross-shore  transects  {x,z)  were 

13  objectively  mapped  with  a  skill  threshold  of  0.95  and  a  signal-to-noise  ratio  of  0.9 

14  (Bretherton  et  al.  1976). 

15 

16  Moorings 

17  Moorings  were  deployed  on  a  4  km  long  cross-shore  transect  at  mean  (tidally- 

18  averaged)  water  depths  H  =S,  15,  18  and  24  m  (M8,  M15,  M18,  M24  in  Eig.  2)  and 

19  instrumented  with  4  to  6  Star-Oddi  (T  only)  or  Sea  Bird  MicroCAT  (S  and  7)  instruments 

20  sampling  at  3  min  intervals.  A  wirewalker,  a  wave-driven  vertically  profiling  platform 

21  (Rainville  and  Pinkel,  2001),  was  deployed  at  7/  =  13  m  (Ml 3  in  Pig.  2),  and 

22  instrumented  with  a  CTD  (Seabird  49,  sample  rate  16  Hz)  and  Chi  a  fluorometer  (WET 

23  Labs  ECO  Triplet,  sample  rate  16  Hz).  Vertical  CTD-i-  Chi  a  M13  profiles  were 

24  completed  approximately  every  2  min;  data  were  averaged  every  0.1  m  in  z  and 
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1  interpolated  onto  regular  4  min  intervals.  Cross-shore  (u)  and  alongshore  (v)  currents 

2  measured  with  acoustic  Doppler  current  profilers  (ADCPs)  deployed  at  the  base  of 

3  moorings  Ml 3  (600  kHz),  M15  (600  kHz)  and  M18  (600  kHz)  were  averaged  at  6  min 

4  intervals  with  1  m  vertical  bins.  A  bottom-mounted  Nortek  Aquadopp  at  M8  (2000  kHz) 

5  sampled  every  2.5  min  with  0.5  m  vertical  bins. 

6 

7  Surf  zone  frames 

8  Fixed  frames  with  acoustic  doppler  velocimeters,  pressure  and  temperature 

9  sensors  (sampled  at  8  Hz,  8  Hz  and  0.2  Hz  respectively)  were  deployed  on  a  cross-shore 

10  transect  spanning  160  m  from  near  the  shoreline  to  4  m  water  depth.  In  addition,  WET 

1 1  Labs  ECO  Triplet  fluorometers  measuring  Chi  a  (sampled  at  0.25  Hz)  were  repeatedly 

12  deployed  for  72  h  periods  on  a  frame  from  within  (Ml. 5,  H  =  1.5  m)  and  seaward  (M4,  H 

13  =  4  m)  of  the  surfzone  (inset  in  Eig.  2b),  facing  30  degrees  from  downward,  nominally 

14  0.5  m  above  the  seafloor.  A  correction  for  data  biases  induced  by  the  turbid  surfzone 

15  conditions  was  applied  the  ECO  Triplets  (Omand  et  al.  2009). 

16 

17  Water  Samples 

18  At  most  casts,  water  samples  (for  Chi  a  analysis  and  phytoplankton  cell  counts 

19  and  identification)  were  collected  near  the  bottom  with  a  messenger-tripped  Niskin  bottle 

20  mounted  1  m  above  the  CTD,  and  at  the  surface  by  hand.  Water  samples  were  also 

21  collected  near  the  surface  every  20  min  on  Oct  12  from  900  to  1500  h  near  Ml. 5  and  M4. 

22  Samples  for  taxonomic  identification  were  preserved  by  transferring  100  ml  of  raw 

23  seawater  to  a  glass  Wheaton  bottle  and  combined  with  10  ml  buffered  37% 

24  Eormaldehyde.  The  preserved  samples  were  analyzed  with  the  Utermohl  settling  method 
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1  (Utermohl,  1958)  and  inspected  at  16X  magnification  under  a  light  microscope.  All 

2  visible  cells  (>5  pm)  were  enumerated  and  identified  to  a  genus  or  species  level  when 

3  possible. 

4 

5  Nonphotochemical  quenching  correction 

6  Nonphotochemical  quenching  (NPQ)  reduces  Chi  a  fluorescence  in  high  light 

7  conditions  (e.g.,  Kiefer  1973,  Falkowski  and  Raven  1997,  Muller  et  al.  2001).  NPQ 

8  corrections,  derived  from  continuous  profiles  of  Chi  a  at  M13  (Fig. 2)  and  surface 

9  irradiance  (measured  with  a  Davis  Vantage  Pro  Plus  cosine  pyranometer),  were  applied 

10  to  all  in  situ  Chi  a  measurements  (see  the  Web  Appendix 

1 1  M13w.aslo.org/lo/toc/vol_xx/issue_x/xxxxal.pdf).  The  NPQ  correction  is  small  and  does 

12  not  affect  the  overall  patterns  and  conclusions. 

13 

14  Results 

15  The  sudden  surface  appearance  of  a  nearshore,  dense  red-tide  patch  is  described, 

16  using  physical  and  biological  data  with  uniquely  high  spatial  and  temporal  resolution. 

17  Nearshore  (H  <  15  m)  Chi  a  was  generally  low  on  5  Oct  (maximum  4.5  pgL'\  Fig.  3a) 

18  and  reached  13  pgL'^  on  8  Oct  and  9  Oct  in  a  subsurface  layer  between  the  16  and  17  °C 

19  isotherms  (Fig.  3b,c).  Maximum  Chi  a  peaked  above  34  pgL'^  on  12  Oct  (Fig.  3d),  and 

20  decreased  to  12  pgL'^  on  16  Oct  (Fig.  3e).  Offshore  (15  m  <  7/  <  50  m)  Chi  a  within  the 

21  subsurface  layer  also  increased,  though  less  intensely,  from  a  Chi  a  maximum  of  6  pgL'^ 

22  on  5  Oct  to  10  pgL'^  on  12  Oct.  Glider  surveys  conducted  approximately  1-10  km 

23  offshore,  concurrently  with  HB06,  also  indicate  an  increase  of  mean  subsurface  Chi  a 

24  between  29  Sept  and  16  Oct  (Todd  et  al.  2009). 
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1  Conditions  on  12  Oct  were  typical  of  a  fall  day  at  Huntington  Beach;  winds  were 

2  weak,  blowing  slightly  offshore  in  the  morning  with  a  northwest  3  ms'^  breeze  in  the 

3  afternoon.  The  surface  water  temperature  was  approximately  17.5  °C,  dropping  to  14  °C 

4  at  the  base  of  the  thermocline  at  ~  20  m  (Fig.  3d).  Offshore  (H  >  15  m  depth),  a  deep  Chi 

5  a  maximum  (DCM)  layer  of  10  pgL'^  Chi  a  lay  between  the  14  and  16  °C  isotherms.  In 

6  the  nearshore  (//  <  15  m  depth)  an  intense  patch  of  Chi  a  was  observed  between  the  16 

7  and  17  °C  isotherms  (Fig.  3d).  Nearshore  and  offshore  phytoplankton  taxonomic 

8  composition  (Table  1)  was  strikingly  different.  The  nearshore  was  dominated  by  a  near 

9  monoculture  of  dinoflagellates  (>  70%  Lingulodinium  polyedrum)  while  offshore  waters 

10  contained  a  greater  diversity  of  genera,  dominated  by  diatoms  from  the  genera 

1 1  Asterionella  and  Chaetoceros.  L.  polyedrum  were  also  the  dominant  nearshore  cells  on  8, 

12  9,  and  16  Oct.  The  nearshore  increase  in  Chi  a  between  5-12  Oct  corresponds  to  a  Chi 

13  a-specific  growth  rate  of  ~  0.3  d'\  This  growth  rate,  reasonable  for  L.  polyedrum  (e.g., 

14  Sullivan  and  Swift,  2003),  suggests  that  the  monotonic  increase  in  subsurface  Chi  a 

15  concentrations  was  largely  generated  by  local  growth  of  this  species. 

16  Five  CTD-i-  Chi  a  cross-shore  transects  and  near-synoptic  jetski  surveys  of  surface 

17  Chi  a  and  T  were  conducted  at  45  min  intervals  between  1200  and  1500  h  in  the 

18  nearshore  region  inshore  of  Ml 3  (Fig.  2).  At  1200  h,  a  high  Chi  a  concentration  (up  to 

19  39  pgL'^)  layer  was  located  ~  5  m  below  the  surface  (Fig.  4a).  Surface  Chi  a  and 

20  temperature  were  spatially  homogeneous  at  3  pgL'^  and  17.7  °C,  respectively  (Fig.  4a,b). 

21  The  hourly-averaged  cross-shore  currents  were  directed  -0.03  ms"'  onshore  near  the 

22  surface  and  offshore  near  the  seafloor.  The  surfzone  alongshore  current  was  0.15  ms'^ 

23  northward.  By  1245  h,  a  small  (-10  m  cross-  and  alongshore)  patch  of  approximately  4 

24  pgL'^  Chi  a  appeared  at  the  surface  with  slightly  cooler,  17.5  °C  water  (Fig.  4c,d),  and  the 
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1  subsurface  high-concentration  Chi  a  patch  moved  inshore  (between  300  <  x  <  600  m). 

2  The  patch  center  lay  just  below  the  17  °C  isotherm. 

3  At  1330  h  (Fig.  4e,f),  a  5  pgL'^  Chi  a  patch  had  further  emerged  at  the  surface 

4  near  x  =  400  m  in  17.5  °C  water,  spanning  100  m  across  shore  and  250  m  alongshore. 

5  The  near-surfzone  region  (x  <  150  m)  had  warmed  to  18  °C.  The  subsurface, 

6  concentrated  Chi  a  patch  was  centered  on  the  17  °C  isotherm  and  increased  in  thickness. 

7  At  1415  h,  the  surface  Chi  a  patch  was  developed,  with  a  7  pgL'^  maximum,  a  cross- 

8  shore  span  of  about  200  m  (Fig.  4g)  and  visually  estimated  alongshore  length  ~  1  km. 

9  The  near-surfzone  region  continued  to  warm  and  the  colder  (<  17  °C)  surface  water  at  x  > 

10  500  m  receded  offshore  (Fig.  4h).  The  Chi  a  patch  continued  to  spread  vertically  into 

1 1  warmer  17.5  to  18  °C  water.  Later  at  1500  h,  the  surface  Chi  a  patch  had  advected 

12  slightly  offshore  with  the  <  17.5  °C  water  (Fig.  4i,j). 

13 

14  Stages  of  the  surface  red  tide  patch  formation 

15  The  sudden  appearance  of  the  red  tide  at  the  surface  was  the  culmination  of 

16  several  physical-biological  processes.  Shoreward  propagation  of  an  internal  wave  trough, 

17  advection  of  a  subsurface  red  tide  layer,  local  formation  of  a  dense  Chi  a  patch  within  the 

1 8  trough,  vertical  patch  spreading,  all  contributed  to  the  appearance  of  the  surface  Chi  a 

19  band.  The  observations  underlying  each  stage,  and  the  subsequent  disappearance  from 

20  very  shallow  water,  are  described. 

21 

22  Shoreward  propagation  of  a  supertidal  internal  wave 

23  On  12  Oct.  between  0000  (midnight)  and  0300  h,  mid- water  column  isotherms  at 

24  M24  moved  upwards  about  5  m,  followed  by  an  approximately  7  m  drop  into  an  internal 
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1  wave  trough  at  0500  h  (star,  Fig.  5a).  The  shoreward-propagating  supertidal  internal 

2  wave  (sIW)  trough  subsequently  passed  Ml 8  and  Ml 5,  reaching  M13  at  1 100  h  (black 

3  stars,  Fig.  5b,c,d).  At  M8  (H  =  7.8  m),  the  sIW  was  obscured  with  no  obvious  isotherm 

4  trough  (Fig.  5e). 

5  The  sIW  propagation  time-lag  At,  the  time  between  arrivals  of  the  sIW  trough 

6  (isotherm  minima,  black  stars.  Fig.  5)  is  estimated  for  adjacent  moorings  (not  including 

7  M8).  These  At  give  the  maximum  correlation  (r  >  0.63,  p  <  0.001,  Table  2)  between 

8  mid-water  column  isotherm  depth  anomalies  (e.g.  Lerczak,  2000)  over  a  15  h  period 

9  (from  0000  to  1500  h)  that  includes  the  sIW  trough  (Fig.  6a).  The  observed  cross-shore 

10  phase  velocity  Cp  =  Ax/At,  where  Ax  is  the  cross-shore  separation  between  adjacent 

11  moorings  decreased  from  0.16  ms'^  between  M24  and  Ml 8  to  0.12  ms"'  between  M15 

12  (black  points.  Fig.  6b).  With  constant  buoyancy  (AO  and  wave  (co)  frequencies,  the 

13  theoretical  phase  speed  (Cp)  for  a  cross-shore  propagating,  mode  1,  internal  wave  (e.g., 

14  Gill,  1982) 

15  (1) 

16  depends  linearly  on  the  depth  H.  Cp  (1)  is  estimated  using  the  vertically-averaged  N 

17  (estimated  from  dT/dz  adjusted  to  density),  and  the  mean  H,  of  adjacent  moorings.  The 

18  buoyancy  frequency  N  (~  0.02  s'^)  was  large  compared  to  co,  so  theoretical  Cp  was 

19  insensitive  to  the  choice  of  co  (Cp  varied  by  <  1%  for  wave  periods  between  2  and  12  hrs: 

20  a  nominal  period  of  2  hrs  was  used).  The  observed  and  theoretical  (1)  Cp  were  similar 

21  (compare  gray  and  black  curves  Fig.  6b).  The  sIW  propagated  coherently  about  3km 

22  onshore  between  M24  and  M13  and  the  observed  Cp  decreased  with  decreasing  depth,  as 

23  often  observed  for  supertidal  internal  waves  on  the  Southern  California  shelf  (e.g., 

24  Winant,  1974;  Winant  and  Olsen,  1976;  Lerczak,  2000;  Johnson  et  al.  2001). 
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Advection  of  the  subsurface  red  tide  layer  and  dense  Chi  a  patch 


2  On  11-12  Oct  at  M13  (H  =  13  m),  the  tide  varied  ±  1  m  (Fig.  7a).  The  cross- 

3  shore  currents  at  Ml 3  were  dominated  by  baroclinic  modes  with  both  diurnal  and 

4  semidiurnal  variability  (Fig.  7b),  consistent  with  prior  studies  at  this  location  (Boehm  et 

5  al.  2002;  Noble  et  al.  2009).  The  subsurface  elevated  Chi  a  (-10-20  pgL'\  Fig.  7c)  layer 

6  repeatedly  shoaled  from  near-bottom  to  near-surface  and  back  down.  The  Chi  a  layer 

7  generally  followed  the  16.0  -  17.0  °C  isotherms  (black  contours  in  Fig.  7c)  throughout  the 

8  diurnal  cycle,  except  for  a  few  hours  after  the  17.5  °C  isotherm  descended,  leaving  a 

9  near-surface  (z  >  -5  m),  5  pgL'^  Chi  a  band  above  the  17.5°  C  isotherm  (e.g.,  at  t  =  -32,  -8 

10  and  14  h  in  Fig.  7c). 

11  At  M13  between  1000  and  1200  h  on  12  Oct,  the  L.  polyedrum  dominated  Chi  a 

12  patch  (Fig.  3d),  maximum  concentrations  was  31  pgL'^  between  5  and  7  m  below  the 

13  surface,  at  the  17.0  °C  isotherm  (Fig.  8).  The  onshore  currents  (approximately  0.05  ms"' 

14  between  1000  and  1200  h)  within  the  layer  (-7  <  z  <  -5  m)  transported  the  patch 

15  shoreward.  Based  on  the  onshore  u  over  the  2  hr  M13  patch  duration,  the  high 

16  concentration  subsurface  Chi  a  patch  was  estimated  to  be  200  m  wide  across-shore, 

17  consistent  with  the  jetski  observations  at  the  surface  (Fig.  6).  Below  the  patch  (z  <  -9  m, 

18  Fig.  8),  Chi  a  fell  to  6  pgL'\  and  the  cross-shore  currents  u  were  directed  offshore.  After 

19  1200  h,  the  upper  layer  Chi  a  concentrations  were  reduced  to  around  10  pgL'^  as  u  in  that 

20  layer  weakened  (<  0.01  ms'^  onshore).  Between  1400  and  1600  h,  the  lower  layer  (z  <  - 

21  7m)  contained  elevated  Chi  a  that  was  transported  offshore  past  M13  by  offshore  u  (Fig. 

22  8).  At  1600  h,  the  upper  layer  (z  >  -7  m)  currents  had  become  offshore  directed, 

23  transporting  a  surface  patch  of  >  4  pgL'^  Chi  a  seaward.  By  1800  h,  Chi  a  had  decreased 

24  to  moderate  (5  pgL'^)  levels  throughout  most  of  the  water  column.  At  2000  h,  a  second 
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1  pulse  of  near-bed  elevated  Chi  a  (near  15  pgL'^)  was  advected  shoreward  as  the  17.0  °C 

2  isotherm  rose  to  z  =  -8  m. 

3 

4  Formation  of  the  dense  subsurface  Chi  a  patch  within  an  internal  wave  trough 

5  Phytoplankton  advection  by  IWs  in  combination  with  depth-keeping  swimming 

6  behavior  may  locally  enhance  Chi  a  above  wave  troughs  as  the  isotherm  depth  beneath 

7  the  surface  Az  increases  (Lennert-Cody  and  Franks  1999,  2002).  Evidence  of 

8  systematically  elevated  Chi  a  within  the  sIW-associated  isotherm  troughs  is  presented, 

9  qualitatively  consistent  with  internal  wave-induced  convergence  and  phytoplankton 

10  depth-keeping. 

11  AtM13onl2  Oct,  Chi  a  in  the  sIW  trough  (black  star.  Fig.  5d),  was  maximum 

12  (31  pgL'^)  and  nearly  triple  the  mean  concentration  (11.5  pgL"')  bordering  the  17.0  °C 

13  isotherm.  Chi  a,  vertically-averaged  between  the  17.0  °C  isotherm  and  the  surface  (solid 

14  black  curve  in  Fig.  5d)  and  the  isotherm  separation  Az  were  significantly  positively 

15  correlated  (Fig.  9a).  Mean  Chi  a  was  maximal  where  the  isotherms  were  the  deepest.  In 

16  contrast,  T,  vertically  averaged  between  the  17.0  °C  isotherm  and  the  surface,  was 

17  uncorrelated  with  Az.  The  linear  relationship  between  the  vertically  averaged  Chi  a  and 

18  Az  and  the  absence  of  a  T-Az  relationship  do  not  depend  on  the  choice  of  isotherm  within 

19  the  range  16.5  to  17.2  °C,  associated  with  the  sIW.  Although  NPQ  can  induce  apparent 

20  along-isotherm  Chi  a  variability  on  a  depth- varying  isotherm,  the  maximum  NPQ 

21  correction  (see  Web  Appendix)  for  Chi  a  on  the  17.0  °C  isotherm  is  9%,  far  smaller  than 

22  the  observed  Chi  a  variability. 

23 

24  Vertical  spreading  of  the  dense  subsurface  Chi  a  patch 
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1  The  surface  Chi  a  patch  appeared  around  1245  h  on  an  isotherm  warmer  {T  =  17.5 

2  °C)  than  the  isotherms  associated  with  the  subsurface  patch  (16  -  17  °C,  Fig.  7c), 

3  possibly  resulted  from  vertical  mixing  of  the  subsurface  Chi  a  patch  with  warmer  surface 

4  water.  Here,  the  role  of  vertical  mixing,  including  cross-isotherm  mixing,  in  the  surface 

5  appearance  Chi  a  patch  is  investigated  by  examining  the  temporal  evolution  of  the  T-z 

6  and  Chi  a-T  relationships  from  1200  to  1500  h,  within  the  750  m-wide  region  between 

7  //  =  5  to  12  m.  This  region  bounds  the  Chi  a  patch  (Fig.  4),  includes  the  estimated  time 

8  and  location  of  sIW  breaking  (white  star  in  Fig  5e),  and  excludes  the  warm  shallow  water 

9  (x  <  200  m,  H  <  5  m.  Fig.  4).  To  minimize  the  effects  of  cross-shore  advection,  CTD-i-Chl 

10  a  cast-derived  Chi  a  and  T  for  each  nearshore  survey  on  12  Oct  were  cross-shore- 

1 1  averaged  (denoted  (Chi  a(z)}  and  (T(z)}),  into  1  m- thick  vertical  bins. 

12  The  depth-average  of  (Chi  a  (z)}  remained  roughly  constant  (within  15%) 

13  between  1200  and  1500  h,  consistent  with  small  depth-integrated  cross-shore  Chi  a 

14  fluxes,  and  small  net  phytoplankton  growth  and  mortality  (not  shown).  The  depth- 

15  averaged  (T(z)}  warmed  slightly,  from  16.9  to  17.1  °C,  perhaps  from  solar  heating.  For 

16  each  survey  (T(z)}  were  approximately  linear  with  depth  (not  shown).  The  vertical 

17  temperature  gradient  d(T(z)}/dz  decreased  monotonically,  and  the  corresponding 

18  buoyancy  frequency  decreased  from  0.018  s'^  at  1200  h  to  0.016  s’'  at  1500  h,  suggesting 

19  that  vertical  spreading  of  T(z),  dominated  the  stratifying  effects  of  solar  warming  over  the 

20  afternoon. 

21  In  addition  to  a  reduction  in  N  (implying  vertical  mixing),  cross-isotherm  Chi  a(z) 

22  mixing  may  also  have  occurred.  For  example,  cross-isotherm  Chi  a  spreading  (and/or 

23  surface  solar  heating)  may  explain  the  appearance  of  the  surface  Chi  a  patch  on  the  17.5 

24  °C  isotherm  at  1330  h  (Fig.  4e,g,i),  since  Chi  a  was  not  present  on  this  isotherm  at  1200  h 
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1  (Fig.  4a).  Subsurface  cross-isotherm  Chi  a  spreading  was  evident  in  the  evolution  of  the 

2  cross-shore-averaged  (Chi  a(z))-(T(z))  relationship  (Fig.  10).  For  all  surveys,  (Chi  a{z)) 

3  was  maximum  near  T  =  16.9  °C.  At  1200  h,  (Chi  aiz))  was  concentrated  over  a  narrow 

4  range  of  T(z),  and  throughout  the  afternoon,  (Chi  a{z))  spread  over  a  larger  range  of  T(z) 

5  (e.g.,  compare  the  dashed  gray  with  the  solid  black  curve  in  Fig.  10).  Although  the 

6  CTD-i-Chl  a  surveys  exclude  the  top  1  m  of  data,  the  increased  (Chi  a{z))  at  {T(z))  >  17.3 

7  °C  after  1200  h,  was  in  qualitative  agreement  with  the  surface  manifestation  of  the  patch 

8  (e.g..  Fig.  4e-j). 

9 

10  Disappearance  of  Chi  a  in  very  shallow  water 

1 1  During  the  jetski-CTD  observation  period,  the  high  Chi  a  surface  patch  did  not 

12  enter  the  surfzone  (0  <  x  <  150  m)  where  the  near-surface  water  in  the  surfzone  warmed 

13  relative  to  near-surface  offshore  waters  (see  Fig.  4).  Between  1000  and  1200h,  near- 

14  surface  T  was  uniform  at  roughly  17  °C  from  within  the  surfzone  to  seaward  of  M13  (Fig. 

15  11a).  Coincident  with  the  Ml 3  subsurface  Chi  a  peak  (Fig.  3d),  M4  Chi  a  (at  0.5  m 

16  above  the  seafloor)  was  elevated  (between  5-10  pgL'^),  with  the  peak  M4  Chi  a  only  ~ 

17  25%  of  that  at  M13  (Fig.  8),  and  varied  on  1  h  time-scales  (solid  black  line.  Fig.  1  lb). 

1 8  The  M4  dinoflagellate  (mainly  L.  polyedrum)  cell  concentrations  also  reached  a 

19  maximum  of  7.3  x  lO"^  cellsL'^  near  1 130  h  (solid  black  line.  Fig.  1  Ic),  while  diatom 

20  abundance  remained  low  0(10  )  cells  L'  (dashed  black  line.  Fig.  1  Ic).  The  bottle  sample 

21  collection  rate  was  not  rapid  enough  to  resolve  sub-hourly  fluctuations.  However,  the 

22  qualitative  agreement  between  dinoflagellate  cell  concentration  and  Chi  a  indicates  that 

23  the  M4  Chi  a  variation  was  due  to  variation  in  total  dinoflagellates.  Within  the  surfzone, 

24  Chi  a  was  lower  at  about  3-4  pgL'^  (Fig.  1  lb),  consistent  with  the  jetski-observed  Chi  a 
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1  values  at  1200  h  (Fig.  4a).  In  the  surfzone,  dinoflagellate  (solid  gray  lines)  and  diatom 

2  (dashed  gray  lines)  concentrations  were  low,  usually  <  lO"^  cells  L‘\  and  relatively 

3  constant  throughout  the  sample  period  (Fig.  1  Ic). 

4  At  1 130  h,  the  M4  Chi  a  began  dropping  and  by  1230  h  reached  the  level  at  Ml. 5 

5  (3  pgL'^).  Ml. 5  and  M4  Chi  a  remained  constant  throughout  the  afternoon  as  T  at  both 

6  locations  rose  (Fig.  1 1).  During  mid-day,  the  near-surface  M13  T  also  increased  slightly 

7  (thin  curve  in  Fig.  1  la).  By  1500  h,  a  cross-shore  temperature  gradient  (with  AT  ~  0.8  C) 

8  had  developed  with  warmest  water  in  the  surfzone  and  cooler  water  offshore,  consistent 

9  with  surface  T  maps  (Fig.  4j).  The  strong  afternoon  T  gradient,  and  morning  gradients  in 

10  Chi  a,  total  cells  and  community  composition  were  maintained  over  only  80  m  in  the 

1 1  cross-shore. 

12 

13  Discussion 

14  A  dense  patch  containing  3X  the  mean  concentration  of  the  red- tide  layer  was 

15  observed  at  M13  (1100  h.  Fig  8).  The  patch  appears  to  have  developed  though  local 

16  intensification  through  the  interaction  of  depth-keeping  swimming,  with  a  shoreward- 

17  propagating  internal  wave  trough.  The  supertidal  internal  wave  arrival  in  shallow  water 

18  (//  <  12  m)  coincided  with  vertical  spreading  of  the  Chi  a  patch,  creating  an  alongshore 

19  Chi  a  band  visible  at  the  surface.  Mechanisms  underlying  these  processes  are  discussed. 

20 
21 
22 

23 

24 
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1  Formation  of  the  dense  subsurface  Chi  a  patch  within  a  shoreward-propagating  wave 

2  trough 

3 

4  A  supertidal  internal  wave  of  depression  (sIW)  propagated  onshore  with  a  phase 

5  speed  consistent  with  linear  theory.  The  sIW  propagated  coherently  past  the  mooring 

6  array  (Fig  5),  maintaining  shape  as  it  moved  into  shallower  waters.  The  wave  passage  at 

7  Ml 3  coincided  with  a  concentrated  patch  of  Chi  a  above  the  sIW  trough  and  within  the 

8  subsurface  Chi  a  maximum  layer.  The  patch  may  have  formed  through  the  combined 

9  effects  of  localized  phytoplankton  depth-keeping  swimming  and  the  IW-induced 

10  convergent  flow  (e.g.  Lennert-Cody  and  Franks, 1999).  Though  propagating  nonlinear 

1 1  internal  waves  may  transport  particles  in  the  direction  of  propagation  above  the  wave 

12  trough  (e.g.  Lamb,  1997),  the  sIW  phase  velocity  (between  0.015  to  0.012  ms'^)  was 

13  significantly  larger  than  the  mean  Eulerian  velocity  within  the  wave  trough  at  each 

14  mooring  (between  0.07  to  0.03  ms'^  at  M24  to  M13  respectively).  The  Chi  a  patch  was 

15  likely  formed  by  local  intensfication  of  the  subsurface  Chi  a  layer,  with  limited  advection 

16  by  relatively  weak  currents,  as  opposed  to  being  advected  from  offshore  within  the  wave 

17  trough.  The  high-concentration  Chi  a  patch  was  estimated  to  be  200  m  wide  across  shore. 

18  This  is  consistent  with  the  cross-shore  width  of  the  surface  Chi  a  patch  (Fig.  4). 

19  The  relationship  between  the  vertically  averaged  Chi  a  and  depth  of  the  17  °C 

20  isotherm  Dz  (Fig.  8a)  is  consistent  with  theoretical  predictions  and  previous  observations 

21  (Lennert-Cody  and  Franks  1999,  2002):  Chi  a  was  enhanced  in  wave  troughs  where 

22  isotherms  are  strained.  The  amount  of  Chi  a  enhancement  depends  on  the  phytoplankton 

23  swimming  proficiency  relative  to  rate  of  the  vertical  movement  of  the  isotherms.  Typical 

24  instantaneous  swimming  speeds  of  L.  polyedrum  in  the  laboratory  range  between  100  and 
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400  mms"'  (Lewis  and  Hallet,  1997;  Buskey,  1997;  Mayali  et  al.  2008).  Vertical  velocities 
(w)  estimated  from  vertical  displacements  of  the  17.0  °C  isotherm  ranged  between  about 
300-1500  mms'^  (c-g-,  2-9  cm  min'\  see  bold  isotherm  displacements  between  1000  and 
1500  h,  Fig.  5d).  The  range  in  w  partially  overlaps  the  range  of  L.  polyedrum  swimming 
speeds,  suggesting  that  the  dinoflagellates  could  have  partially  or  completely  countered 
the  sIW-generated  vertical  velocities. 

Local  Chi  a  enhancement  above  isotherm  troughs  indicates  that  Chi  a  is  not  a 
passive  tracer  in  this  region  of  large  amplitude,  shoreward-propagating,  internal  waves. 
Swimming  may  contribute  significantly  to  deviations  of  the  Chi  a  concentration  from 
predictions  of  a  passive  tracer,  but  quantification  requires  more  detailed  understanding  of 
both  L.  polyedrum  swimming  behavior,  and  of  potential  Lagrangian  transport  in  nonlinear 
internal  waves. 


sIW  breaking  in  shallow  water 

The  sIW  trough  was  not  evident  at  M8  (H  =  8  m,  white  star.  Fig.  5e),  possibly  due 
to  IW  breaking  (e.g.,  Cacchione  and  Wunsch,  1974).  At  M13  (H  =  13  m),  the  sIW 
amplitude  ao  was  ~  5  m,  based  on  the  vertical  displacement  (crest  to  trough)  of  the  17  °C 
isotherm  (thick  black  contour.  Fig.  5d).  The  group  velocity  of  a  shoreward-propagating 
internal  wave  decreases,  and  ao  increases,  to  conserve  wave  energy,  until  the  wave  breaks 
(e.g..  Gill,  1982).  For  an  idealized  2-layer  system,  breaking  occurs  when  the  undisturbed 
lower  layer  thickness  Diow  is  less  than  Dbp  (Helfrich,  1992)  defined  as 

a 

D  = - ^ -  (2) 

(1.6  i  + 0.6)’ 


19 


1 


where  A  is  the  wavelength  and  L  is  the  distance  to  the  shelf  break;  L  =  8000  m  at 


2  Huntington  Beach,  where  //  =  60  m  depth.  The  ratio  A/L  is  small  for  supertidal  IWs,  and 

3  Dbp  depends  only  weakly  on  A.  Using  the  observed  Cp  =  0.10  ms’'  and  A  ~  500  m)  (given 

4  a  wave  period  ~  2  hrs)  an  idealized,  sIW  with  ao  =  5  m,  would  begin  breaking  when  Di„„ 

5  <  Dbp  ~  7  m. 

6  Taking  the  lower  layer  thickness  Diow  as  the  elevation  of  the  17.0  °C  isotherm 

7  above  the  seafloor,  Diow  was  11m  before,  and  8  m  after,  the  trough  passed  Ml 3  (Fig.  5d) 

8  and  the  sIW  was  not  breaking  according  to  (2).  At  M8  (H  =  7.8  m),  Diow  was  always  less 

9  than  Dbp  (Fig.  5d,e)  indicating  that  wave  breaking  has  occurred  according  to  (2). 

10  Helfrich's  criteria  (2),  and  the  observed  obscuring  of  the  sIW  trough,  indicates  that  IW 

1 1  breaking  occurred  between  M13  and  M8.  A  projection  of  the  linear  relationship  between 

12  Cp  and  H  from  (1)  (Fig.  6b),  predicts  that  the  sIW  likely  began  breaking  near  1230  h 

13  (white  star  in  Fig  5e). 

14 

15  Vertical  mixing  of  Chi  a 

16  On  12  Oct,  the  peak  horizontally-integrated  nearshore  Chi  a  decreased  over  time 

17  and  Chi  a  was  spread  over  a  broadening  range  of  isotherms  (Fig.  10).  The  shoreward- 

1 8  propagating  IW  breaking  was  anticipated  in  ~  8  m,  and  may  have  mixed  the  water 

19  column  vigorously  enough  for  the  subsurface  bloom  to  reach  the  surface,  where  it  was 

20  observed  visually  (Fig.  1)  and  with  the  jetski  maps  (Fig.  4c,e)  as  a  200  m-wide,  surface 

21  patch  onshore  of  M13.  The  vertical  spreading  of  the  subsurface  (Chi  a{z))  layer,  which 

22  includes  cross-isotherm  mixing,  is  quantified.  The  (Chi  a{z))  center  of  mass  ((z)chia)  (e.g., 

23  Csanady,  1973), 
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deepened  linearly  from  -7.1  m  to  -7.8  m  (Fig.  12a)  over  3  hrs  at  a  rate  of  60 
pms'V  The  deepening  of  the  layer  peak  is  qualitatively  evident  in  the  CTD-i-  Chi  a 
sections  (Fig.  4a,c,e,g,i).  The  phytoplankton  layer  squared  half-width  (o^ chiaXe-.g., 
Csanady,  1973), 


J(z-{z}chjHChla(z))dz 

C/  Chla  =  — - 5 - ,  (4) 

jf  {Chla(z))dz 

-H 


increased  roughly  linearly  with  time  until  slowing  around  1500  h  (solid  line,  Fig.  12b),  as 
anticipated  with  breaking  sIW  turbulence.  The  layer  width  (2o)  varied  between  4.2  m  at 
1200  h  and  4.7  m  at  1415  h,  consistent  with  the  qualitative  Chi  a  vertical  spreading  seen 
in  the  CTD-i-Chl  a  sections  (Fig.  4a,c,e,g,i).  The  slowing  of  the  Chi  a  vertical  spreading 
by  1500  h  may  be  attributed  to  subsiding  of  the  breaking  sIW-generated  mixing.  The 
vertical  diffusivity  Kzz  and  squared  layer  width  c/  are  related  (e.g.,  Csanady,  1973), 


do‘ 


K 


Chla 


zz  2 


dt 


(5) 


The  best-fit  linear  slope  of  chia  between  1200  and  1415  h  (dashed  line,  Fig.  12b)  yields 
Kzz  =  0.5  X  10'"^  m^s'\  close  to  the  range  for  diffusion  of  surface-layer  (nitrate  and 
temperature)  of  0.6  x  10'"^  to  3  x  lO'"*  m^s'^  previously  estimated  on  the  Southern 
California  shelf  (e.g.,  Eppley  et  al.  1979;  Dillon  and  Caldwell,  1980). 

Phytoplankton  swimming  behavior  is  implicated  in  IW  trough  Chi  a 
intensification  at  Ml 3  (Fig.  9a).  The  60  pms'^  downward  velocity  of  the  Chi  a  patch 
center  of  mass  (z)chia  between  1200  and  1415  h  was  small  compared  with  typical  L. 
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1  polyedrum  swimming/sinking  speeds  (100  -  400  pms'^),  and  swimming  or  advection 

2  could  have  caused  the  downward  movement  of  (z)chia  between  1200  and  1500  h.  The 

3  potential  Chi  a  vertical  spreading  rate  owing  to  random  phytoplankton  swimming  (akin  to 

4  a  random  walk)  can  be  scaled  by  where  Us  is  a  typical  swimming  speed,  and  x  the 

5  decorrelation  time  of  swimming  direction.  Assuming  Us  ~  200  pms‘\  t=  40  min  is 

6  required  to  reproduce  the  observed  Kzz.  This  x  is  implausible  large,  and  would  induce 

7  ballistic  dispersion  (ie.  chia  ~  Csanady,  1973),  instead  of  the  observed  linear 

8  (Brownian)  relationship  (Fig.  12b).  The  vertical  diffusion  of  the  Chi  a  layer  is  not  due  to 

9  random  phytoplankton  swimming,  and  was  likely  induced  by  breaking  internal  wave 

10  turbulence.  The  expected  sIW  arrival  time  (-1230  h)  based  on  (1),  and  the  anticipated 

1 1  breaking  of  the  IW  onshore  of  Ml 3  (based  on  (2),  together  with  the  evidence  for  vertical 

12  mixing  of  T  and  Chi  a,  suggest  that  sIW  breaking  drove  the  appearance  of  the  surface 

13  alongshore  Chi  a  band. 

14 

15  A  warm-water  barrier  to  the  Chi  a  patch  in  the  surf  zone 

16  The  surface  Chi  a  patch  did  not  enter  the  surfzone.  A  strong  gradient  in  T,  Chi  a 

17  and  phytoplankton  was  observed  over  80  m  spanning  the  surfzone  (Fig.  1 1).  The  warm  > 

18  18  °C  surfzone  water  present  at  Ml. 5  and  M4  after  1300  h  likely  originated  from  the 

19  Talbert  Marsh,  located  1  km  South  (+y)  of  the  fixed  instrument  array  (see  Fig.  2), 

20  creating  a  barrier  to  Chi  a  entry  into  the  surfzone  from  offshore.  Before  1200  h,  the  flood 

21  tide  (Fig.  7a)  raised  the  marsh  water  level.  Over  the  morning  and  early  afternoon,  the 

22  water  within  the  shallow  marsh  (//  -  2  m)  was  solar  heated  relative  to  near-surface 

c  o 

23  temperatures  at  M13.  Using  the  marsh  area  (10  m  ,  Grant  et  al.  2001),  tidal  prism  (2.35  x 

24  10  m,  Jeong  et  al.  2008)  and  the  solar  insolation,  the  marsh  water  would  be  heated 
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1  about  0.8  °C.  With  the  ebb  tide  at  1200  h,  the  warmed  marsh  water  began  draining  into 

2  the  surfzone.  The  0.15  ms'^  northward  surfzone  alongshore  current  (see  Fig.  4),  driven  by 

3  obliquely  incident  breaking  surface  gravity  waves  (e.g.,  Longuet- Higgins,  1970), 

4  transported  this  water  to  the  study  region.  The  predicted  increase  of  -0.8  °C  at  Ml. 5  and 

5  M4  is  roughly  consistent  with  the  near-surface  T  gradient  that  developed  between  the 

6  surfzone  Ml. 5  and  M13.  The  T increase  coincided  with  the  reduction  of  Chi  a  at  M4 

7  (Fig.  11)  and  suppressed  exchange  between  the  surfzone  and  the  water  seaward.  If  the 

8  warm  water  had  not  been  present,  the  intense  Chi  a  patch  may  have  entered  the  surfzone. 

9 

10  Summary  and  synthesis 

1 1  Observations  and  mechanisms  underlying  the  sudden  surface  appearance  of  a 

12  nearshore  red  tide  have  been  described.  The  kilometer-long  shore-parallel  surface  Chi  a 

13  band  (maximum  Chi  a  =  1  pgL'^)  persisted  for  approximately  3  hrs  at  //  =  6-8  m  with  a 

14  cross-shore  width  of  200  m.  A  subsurface  red  tide  layer  (dominated  by  Lingulodinium 

15  polyedrum)  developed  in  the  nearshore  over  4-7  days  prior  to  the  appearance  of  the  patch 

16  at  the  surface.  At  the  M13,  an  elevated  Chi  a  patch  was  observed  in  the  trough  of  a 

17  supertidal  internal  wave  (sIW)  that  propagated  coherently  across  the  mooring  array  (Fig. 

18  13a).  The  elevated  Chi  a  in  the  sIW  trough  is  consistent  with  phytoplankton  swimming 

19  and  IW-induced  convergent  flow  intensifying  Chi  a.  Based  on  the  estimated  phase  speed, 

20  the  sIW  arrived  in  shallow  water  (<  12  m)  around  1230  h,  where  wave  breaking  was 

21  anticpated,  generating  turbulence  that  led  to  vertical  mixing  and  the  appearance  of  the 

22  Chi  a  patch  at  the  surface  (Fig.  13b).  The  wave  arrival  time  (1230  h),  the  approximate 

23  patch  width  (200  m),  and  depth  of  anticipated  wave  breaking  (H  -  8  m)  are  consistent 

24  with  the  appearance  and  dimensions  of  the  surface  Chi  a  band  (Fig.  4). 
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1  Between  1400  and  1700  h,  the  lower  portion  of  the  vertically  spread  Chi  a  patch 

2  appears  to  have  been  advected  offshore  past  Ml 3  with  the  baroclinic  cross-shore  currents 

3  (Fig.  8).  The  surface  Chi  a  patch  persisted  at  a  warmer  isotherm  (>  17.5  °C)  than 

4  originally  associated  with  the  Chi  a  patch  (Fig.  13c),  consistent  with  cross-shore  mixing. 

5  The  surface  Chi  a  patch  was  transported  offshore  past  M13  after  1600  h  by  the  offshore- 

6  directed  surface  currents  once  the  baroclinic  internal  tide  had  reversed  (Fig.  13d). 

7  The  drivers  of  Chi  a  distributions  in  the  nearshore  are  complex.  The  interaction 

8  between  motile  phytoplankton  and  tidal  currents,  propagating  internal  waves,  breaking 

9  internal  waves,  and  surface  wave-driven  surfzone  currents  can  create  phytoplankton 

10  patchiness  and  dramatic  spatial  and  temporal  changes  in  community  structure  over  lO’s 

11  of  meters  and  10’ s  of  minutes.  These  mechanisms  that  concentrate,  advect,  and  mix 

12  nearshore  phytoplankton,  changing  exposure  durations  and  intensity,  may  have  large 

13  consequences  for  the  ecology  of  the  nearshore  region. 
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Web  Appendix:  Corrections  for  Non-Photochemical  Quenching  (NPQ) 
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Chi  a  concentration  is  linearly  related  to  Chi  a  fluorescence  under  most  oceanic 
conditions  (Lorenzen,  1966).  However,  in  situ  Chi  a  fluorescenee  measurements  can  be 
biased  low  due  to  a  phytoplanktonic  physiological  adaptation  called  nonphotochemical 
quenching  (NPQ:  Kiefer  1973,  Falkowski  and  Raven  1997,  Muller  et  al.  2001).  NPQ 
protects  cells  from  high  light  conditions  near  the  surface.  To  prevent  cell  damage  from 
light  absorption  beyond  the  saturated  capacity  for  photosynthesis,  a  fraction  of  the  excess 
light  energy  is  converted  by  the  cell  to  heat,  rather  than  fluorescence.  The  result  is 
reduced  fluorescence  capacity  per  cell,  leading  to  underestimation  of  fluorescence- 
derived  near-surface  daylight  Chi  a  concentration. 

Here,  NPQ  correetion  of  fluorescence-derived  Chi  a  was  based  upon  the  cell’s 
estimated  local  light  environment  (following  Hodges,  2006).  A  similar  approach 
presented  in  Todd  et  al.  (2009)  with  a  correction  based  on  local  irradiance,  produced  a 
nearly  identical  result  (maximum  6.9%  difference)  to  the  Chi  a  correction  presented 
below.  The  'corrected'  Chi  a  (z,t)  (in  the  absence  of  NPQ)  is  a  function  of  the 
fluorescence  measured  Chi  a,  Chi  am(z,t),  and  a  quenching  function  Q(z,t),  i.e., 

Chla(z,t)  =  ^^^^f^  ,  (Al) 

Q{z,t) 


where  Q(z,t)  is  a  function  of  local  irradiance  (j)(z,t).  Hodges  (2006)  suggests  a  functional 
form  for  Q(z,t), 


Q(z,t)  = 


k 

k  +  (l)*(z,t)’ 


(A2) 


where  k  is  a  constant,  and  (p*(z,t)  is  the  normalized  local  irradiance.  The  downward 
surface  (z  =  0)  irradiance  0o(t)  was  measured  for  24  hrs  near  the  field  site  on  16  Oct  2006, 
4  days  after  the  red  tide  observations  (12  Oet).  Both  12  Oct  and  16  Oct  were  sunny. 
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1  cloudless  days,  thus  16  Oct  (j)o(t)  was  assumed  representative  of  surface  irradiance 

2  conditions  during  similar  days  at  HB06.  Surface  irradiance  was  normalized  by  the 

3  maximum  surface  irradiance,  i.e.,  =  <po(t)/<pomax-  Irradiance  was  assumed  to  decrease 

4  exponentially  with  depth  according  to  (l)(z,t)  =  (poi t)e^'  where  y  is  the  unknown  decay 

5  scale  estimated  from  the  vertical  length  scale  of  the  covariation  of  Chi  am(z,0  with  ^o*( t). 

6  The  red  tide  time  (5  Oct  through  16  Oct)  period  was  not  used  due  to  the  intense 

7  subsurface  Chi  a  layer  with  large  vertical  excursions  over  a  24  hr  period  (see  Fig.  7c)  that 

8  would  complicate  the  estimation  of  y.  Instead,  the  time  period  from  20-25  Sept  2006, 

9  when  M13  Chi  amizj)  was  low  (1  -  4  pgL'^)  with  no  subsurface  maxima,  was  used  for  the 

10  covariation  analysis.  Chi  amizj)  covaried  with  0o*( 0  at  the  surface  (r^  =  0.50,  p  <  0.01). 

1 1  The  correlation  diminished  with  depth  over  an  e-folding  scale  y  =  9.7  m.  This  y  is 

12  comparable  to  prior  estimates  (y  ~  5  -  1 1  m)  from  coastal  waters  in  the  Southern 

13  California  Bight  (e.g..  Holmes,  1970,  Conversi  and  McGowan,  1994). 

14  In  (A2),  the  optimal  k=  1.6  was  chosen  to  minimize  the  correlation  between  Chi 

15  a(z,t)  and  Q(z,t),  resulting  in  Q(z,t)  ranging  between  0.615  and  1.  Under  maximum 

16  irradiance  [i.e.,  0o*(12OO  h)  =  1],  near-surface  Chi  a  was  enhanced  by  56%  relative  to  Chi 

17  am,  similar  to  Hodges  (2006)  finding  that  the  maximum  near-surface  correction  was 

18  roughly  100%  Chi  am-  The  NPQ  correction  derived  in  (Al)  and  (A2)  was  applied  to  all  in 

19  situ  Chi  am  measurements  (CTD-i-Chl  a,  M13,  Ml. 5,  M4,  see  Fig.  1)  presented  here. 

20  Although  the  NPQ  correction  modified  Chi  a  concentrations  as  much  as  2X  near  the 

21  surface,  the  patterns  described  here  were  also  present  in  uncorrected  data.  For  example, 

22  along  the  17  °C  isotherm  the  maximum  NPQ  correction  imposed  a  9%  change  from  the 

23  raw  measurement;  insufficient  to  explain  the  3X  Chi  a  increase  observed  in  the  sIW 

24  trough. 
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Tables 


Table  1.  Composition  of  phytoplanktonic  genera  on  12  Oct  within  the  nearshore  (H  <  15 
m  depth)  and  offshore  (H>  15  m)  regions.  Dinoflagellates  are  identified  in  bold  type. 


Genus 

%  Nearshore 

%  Offshore 

Lingulodinium 

71.6 

4.9 

Chaetoceros 

6.4 

34.8 

Asterionella 

6.3 

32.9 

Prorocentrum 

4.6 

0.9 

Psuedo-nitzschia 

4.3 

7.1 

Skeletonema 

0.6 

7.5 

Other 

6.1 

11.9 
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-y 

Table  2.  Separation  distance  (Ax),  time  lag  maximizing  r  (At)  and  correlation  coefficient 
(r  )  of  the  IW  depth  contours  (Fig.  9b)  between  neighboring  moorings.  All  correlations 
are  significant  at  the  99%  level. 


Location 

Ax 

(m) 

At  (h) 

M24toM18 

1643 

2.9 

0.75 

M18  toM15 

715 

1.6 

0.63 

M15  toM13 

549 

1.3 

0.79 
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Figures 


/ 


Figure  1 .  Un-retouched  photograph  of  a  reddish  alongshore-parallel  band  of  red  tide 
(arrow)  approximately  500  m  offshore,  taken  from  the  beach  in  early  afternoon  on  12  Oct 
The  small  breaking  surface  wave  in  the  foreground  is  near  the  beach,  and  the  5  m  long 
CTD+Chl  a  cast  boat  is  offshore. 
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Figure  2.  Schematic  of  HB06  instrumentation:  (a)  plan  view  of  bathymetry  contours 
versus  the  alongshore  (y)  and  cross-shore  (x)  coordinates  and  (b)  cross-shore  transect  of 
bathymetry  (dark  shading).  The  vertical  coordinate  is  z,  with  z  =  0  m  at  the  tidally 
averaged  surface  and  positive  upward.  Moored  temperature  strings  were  located  at  total 
depths  H  =8,  15,  18  and  24  m  (circles,  M8  to  M24).  A  vertically  profiling  CTD-i-Chl  a 
wirewalker  (Ml 3)  was  located  at  //  =13  m.  Bottom-mounted  ADCPs  were  located  near 
moorings  M8,  M13,  M15  and  M18  (stars).  Fixed  frames  within  the  surfzone  (Ml. 5)  at  H 
=  1 .5  m  and  seaward  of  the  surfzone  (M4)  at  H  =  4  m  measured  Chi  a  (triangles),  T 
(circles),  wave  height,  and  currents  (see  inset  in  (b)).  A  GPS-tracked  jetski  measured 
surface  T  and  Chi  a  within  the  black  dashed  box  in  (a)  and  small  boat  (CTD-i-Chl  a) 
completed  transects  during  daylight  hours  from  //  =  5  to  60  m  (~8  km  offshore).  The 
Talbert  Marsh  outlet  is  located  ~1  km  south  of  the  instrument  transect. 
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Figure  3.  Cross-shore  sections  (x,z)  of  Chi  a  (|agL'\  colors)  and  T  (°C,  black  contours) 
conducted  roughly  between  0900  and  1 130  h  from  5  Oct  to  16  Oct  .The  red  circle  is  Ml 3 
surface  location.  The  white  line  traces  the  4  |agL''  Chi  a  contour.  The  black  box  indicates 
the  nearshore  red  tide  focus  area  on  12  Oct. 
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Figure  4.  Sequential  (top  to  bottom)  maps  of  nearshore  Chi  a  (left)  and  T  (right)  on  12 
Oct  from  jetski-based  horizontal  {x,y)  measurements  and  CTD+Chl  a  based  cross-shore 
transect  ix,z)  measurements.  Arrows  represent  the  directions  of  the  near-surface 
alongshore  (at  Ml. 5,  M8  and  M13)  and  near-surface  (~2  m  below  mean  sea  level)  and 
near-bottom  (~  2  m  above  the  seafloor)  cross-shore  currents  (at  M8  and  M13  only).  The 
jetski  required  30  min  to  complete  the  8  cross-shore  transects  used  for  each  map,  and  the 
centered-time  (1200-1500  h)  is  shown  to  the  right  of  each  panel.  Each  CTD-i-Chl  a 
transect  consisted  of  5  or  6  casts.  Black  lines  are  temperature  contours.  The  white  line 
traces  the  4  pgL'^  Chi  a  contour. 
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Figure  5.  Contours  of  T  (°C)  versus  time  and  z  at  (a)  M24,  (b)  M18,  (c)  M15,  (d)  M13,  (e) 
M8.  Black  stars  represent  an  internal  wave  trough  (see  bold  contours)  which  can  be  seen 
to  propagate  onshore.  The  white  star  is  the  trough  location  at  M8  predicted  by  linear 
theory.  In  (d)  Ml 3,  Chi  a  (pgL'^)  is  colored.  Time  zero  corresponds  to  0000  h  12  Oct. 
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Figure  6.  (a)  Depth  anomaly  of  mid  water-column  isotherms  (bold  contours  in  Fig.  5) 
versus  lagged  time  {t-At)  at  M24,  M18,  M15,  and  M13  (see  legend),  with  At  for  each 
mooring  chosen  to  place  the  isotherm  trough  minimum  at  t=0.  (b)  Observed  (solid  black 
line)  and  theoretical  phase  velocity  Cp  (solid  gray  line,  eq.  1)  versus  the  water  depth  H. 
The  observed  Cp  =Ax/At  was  estimated  midway  between  neighboring  moorings,  with  Ax 
the  mooring  separation.  The  theoretical  Cp  is  for  a  linear,  mode-one  IW  with  2  h  period, 
and  using  the  depth- averaged  buoyancy  frequency  N  (solid  gray)  plus/minus  standard 
deviation  of  N  (dashed  gray)  between  moorings. 
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Figure  7.  Time  series  (spanning  54  hrs)  at  Ml 3  (H=  13m)  of  (a)  tidal  excursion  from 
mean  sea  level,  (b)  vertical  profiles  of  hourly  averaged  cross-shore  currents  m  at  M13 
(positive  =  onshore)  and  (c)  15  min  averaged  vertical  profiles  of  Chi  a  (pgL'\  colors)  and 
r  (  C,  black  contours)  at  M13.  The  white  curves  trace  the  4  pgL''  Chi  a  contour.  Time 
zero  is  0000  h  12  Oct. 
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Figure  8.  Time  series  (spanning  11  hrs  of  Fig.  7)  of  15-min  averaged  vertical  profiles  of 
Chi  a  (pgL'\  colors)  and  T(C,  black  contours)  at  M13  (//  =  13  m).  Time  10  hours 
corresponds  to  1000  h  on  12  Oct.  The  hourly  averaged  cross-shore  currents  (u)  are 
represented  by  arrows  for  \u\  >  0.01  ms"'  and  by  black  dots  where  ImI  <  0.01  ms"'.  The 
white  line  is  the  4  pgL'^  contour. 
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Figure  9.  Vertically  averaged  (a)  Chi  a  and  (b)  T  above  the  17.0  °C  isotherm  versus 
isotherm  depth  {Az)  at  4  min  intervals  between  0600  and  1700  h  (number  of  data  points,  n 
=  165)  on  12  Oct  2006.  In  (a),  the  best  linear  fit  (dashed  line)  has  slope  1.29  ±  0.03  pgL' 
*m'',  and  /  =  0.57,  p  <  0.01.  In  (b)  /  =  0.013,  p  >0.1. 
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Figure  10.  Chi  a  versus  T  (at  1200,  1330,  1415,  and  1500  h)  based  upon  the  CTD+Chl  a 
surveys  (between  //  =  5  to  12m  depth,  see  Fig.  4)  horizontally  averaged  along  1  m 
isobars. 
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Figure  11.  Time  series  from  0930  to  1800  h  on  12  Oct  of  (a)  T  (1  min  averages),  (b)  Chi  a 
(1  min  averages)  and  (c)  dinoflagellate  (solid)  and  diatom  (dashed)  concentrations  from 
M4  {H=  4m)  and  Ml  .5  (H  =  2  m),  and  in  (a)  20  min  average  T,2m  below  the  mean  sea 
surface  at  M13  (H  =  13  m)  (see  legend).  The  CTD+jetski  sampling  period  (1200  to  1500 
h,  Fig.  4)  is  indicated  by  the  gray  bar  between  vertical  dashed  lines. 
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Figure  12.  Time  series  from  1 130  to  1530  h  on  12  Oct  based  upon  the  CTD+Chl  a 
surveys  (between  //  =  5  to  12  m  depth,  see  Fig.  4)  of  (a)  depth  {z)chia  of  the  center  of 
mass  of  cross-shore  averaged  (Chi  a  {z))  and  (b)  squared  half-width  c/ chia  (eq.  4)  of  the 
Chi  a{z)  layer.  In  (a),  the  linear  best  fit  (dashed  line)  yields  a  downward  velocity  of  60 
pms'^  and  in  (b)  the  linear  best  fit  (over  the  first  four  data  points,  dashed  line)  gives  a  Chi 
a{z)  vertical  diffusivity  ~  0.5  x  lO'"^  m^s’'. 
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Figure  13.  Schematic  summarizing  the  12  Oct  2006  observations.  The  locations  of  Ml 3 
(red  circle  w/  vertical  line,  1  km  from  shore,  //  =  13  m)  and  the  Ml. 5  and  M4  frames  are 
indicated  (see  triangles  with  vertical  lines.  Fig.  2).  Jetski  and  CTD+Chl  a  casts  provided 
additional  observations,  (a)  Chi  a  in  an  sIW  trough  was  locally  concentrated  by 
swimming  and  sIW  strain  as  the  wave  propagated  onshore,  (b)  sIW  breaking  in  shallow 
water  caused  vertical  mixing  and  outcropping  of  the  17.5  °C  isotherm  and  upper  portion 
of  the  Chi  a  patch,  (c)  deeper  Chi  a  was  advected  offshore  with  the  near-bottom 
baroclinic  currents  and  (d)  surface  Chi  a  is  advected  offshore  after  the  reversal  of  the 
baroclinic  currents. 
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